Master Equation and Control of an Open Quantum System with Leakage 
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Given a multilevel system coupled to a bath, we use a Feshbach P, Q partitioning technique to 
derive an exact trace-nonpreserving master equation for a subspace Si of the system. The resultant 
equation properly treats the leakage effect from Si into the remainder of the system space. Focusing 
on a second-order approximation, we show that a one-dimensional master equation is sufficient to 
study problems of quantum state storage and is a good approximation, or exact, for several analytical 
models. It allows a natural definition of a leakage function and its control, and provides a general 
approach to study and control decoherence and leakage. Numerical calculations on an harmonic 
oscillator coupled to a room temperature harmonic bath show that the leakage can be suppressed 
by the pulse control technique without requiring ideal pulses. 
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Introduction. — Control of quantum dynamics is of 
great interest for "quantum technology industries" , such 
as quantum computing. The control of closed quan- 
tum systems is well established, and has been extensively 
studied in areas such as chemical physics P, 0] ■ Efforts 
to extend these studies to open systems, where the sys- 
tem interacts with an environment, are now underway[3[. 
Quantum information processing has already expended 
considerable effort on open systems. Hence, we antic- 
ipate that methods being developed in the latter area 
may well be useful in the former, and vice- versa 

A fundamentally difficult problem in quantum infor- 
mation processing is that of decoherence, i.e., the loss of 
quantum information in a system due to its interaction 
with its environment (or "bath" ) [1, 0| • For multi- level 
systems, such as molecules, the interaction can also cause 
leakage (i.e. loss of population) from a system subspace 
of interest, denoted Si, into the system space outside 
of Si. Theoretical strategies for combating such dele- 
terious environmental effects in the absence of natural 
decoherence free subspaces d, 0] , invoke the dynamical 
control of system-environment interactions by external 

fields @, S 0, [nl, E [il [11 . 

The aim of dynamical control in open systems is to 
suppress effects of the environment in order to control 
system processes at will. For example, ideal Bang- Bang 



equation for the density operator of the system up to sec- 
ond order in the system bath coupling. This equation al- 
lows the tractable treatment of complicated quantum dy- 
namical processes, eliminating the ideal zero width pulses 
and Trotter formula assumptions and allowing natural 
dynamical evolution and dynamical control on an equal 
footing. For example, Ref. [ls|| provides a unified way 
to suppress the decoherence of two level systems by ar- 
bitrary fields that control the system-bath interaction. 

Below, we consider an N dimensional system {N can 
be infinite), spanned by the bases coupled to a 

bath, and develop procedures to protect quantum infor- 
mation stored in a d-dimensional subspace Si of the sys- 
tem. Specifically, we use a projection operator approach 
to obtain an exact trace-nonpreserving master equation 
for the dynamics of the system subspace Si . We then in- 
troduce a characteristic leakage function which is used to 
consider dynamics and control of this subspace to second 
order in the system-bath interaction. 

Zwanzig's Projection- Operator Approach. — The most 
general Hamiltonian of the N dimensional system plus 
the bath is 



H — Hq + Hj — Hs + Hb + Hi, 



(1) 
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(BB) control of decoherence, decay and leakage 
utilizes idealized zero-width pulses and the Trotter for- 
mula to achieve this goal, although higher-order finite- 
pulse widths have been considered in model cases [l^ . 

Proposals to control decoherence by realistic (non- 
ideal) pulses have often invoked Zwanzig's projection op- 
erator techniques l^, 14 1 resulting in a differential master 



where Hs and Hb are the system and bath Hamiltonians, 
respectively, and the system-bath interaction is Hi — 
X^a SaBa, whcrc Sa and Ba are Hermitian system and 
bath operators. During the dynamics, the information 
stored in the system subspace of interest Si will distribute 
into the bath and leak into the system states outside of 
Si if one does not protect the subspace. Therefore, in 
order to store the information in Si, we need to control 
the system through an interaction with an external field 
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Hc{t), protecting the information during the course of 
time t. The total system Hamiltonian including control 
is therefore Hsit) ^ Hs + H^{t). 

We first derive a closed master equation for Si. Since 
Si opens both to bath and to other system states outside 
of Si , effects of leakage have to be considered in deriving 
a master equation. As usual, the derivation begins in the 
interaction representation with respect to Hq, in which 
the equation of motion is §i p{t) = —i[Hi{t),p{ty\ = 
C{t)p{t) . Here the system-bath interaction is Hi{t) in 
the interaction representation and the LiouviUe superop- 
erator C{t) is defined by this equation [l^. The super- 
projection operation V that we seek defines the relevant 
part of the total density matrix "Ppit) for our new 
open system, i.e. the subspace Si of the entire system. 
Specifically, the superprojection operation comprises two 
commuting parts: a trace over the bath components, 
and a projecting out of the Si subspace from the full 
A^-dimensional space of the system. The associated su- 
perprojector V is therefore defined as 

Vp = P trsipjP ® pB = ri^ pB, 

where P = 'Yl'i=o I*) (*l denotes a projection onto Si and 
•q is the relevant part of the total density matrix, which 
is projected from the total system density matrix ps — 
tiBP as rj = PpsP. The matrix is chosen as the 
initial state of the bath. The trace of the matrix rj 
is not necessarily one, but rj satisfies a closed equation 
as does ps, and plays the same rule as ps, but in Si. 
That is, an arbitrary system observable acting only on 
Si obeys the relation O = POP. The expectation value 
of the operator O when the system-|-bath is in state p is 
tr{pO} =tr{Pp5P0} =tr{?70}, which is the same as the 
expectation value of the operator O in the state charac- 
terized only by rj. This implies that the matrix rj provides 
a complete description of the physics in Si, in the same 
sense that ps completely describes the total open system. 

Applying the superprojector to the equation of motion 
for p gives a time-local master equation 



di 



vp{t) = mvpit), 



(2) 



where JC{t) is the time-convolutionless generator [iq . 
Unlike the usual approach, our new system is a 
d— dimensional subspace of the A'^-dimensional space 
of the total system. Alternatively, Eq. ([2]) can be 
derived by applying the Feshbach projection operator 
approach 1^ 17 1 to the traditional trace-preserving mas- 
ter equation for ps. 

Equation ([2]) is exact and holds for almost [l^ all arbi- 
trary systems and interactions, and for initial conditions 
p{0) = Vp{0), i.e. where the quantum system is initially 
within Si. Since population can flow out of this subspace, 
the master equation is not trace-preserving. Unfortu- 
nately this equation is as difficult to solve as the original 



equation. Therefore, perturbation expansions are needed 
in order to apply the result to actual problems. 

To second order in the coupling strength of the interac- 
tion, JC{t) = /p dsVC{t)£{s)V. Introducing the explicit 
expressions for the projection operator and the LiouviUe 
superoperator, we can obtain the second-order x d di- 
mensional master equation in the interaction representa- 
tion, 

^ Vit) = -A' J dsP trB[Hi{t), [Hi{s), 7^{t) ® pb]P. 

(3) 

Here \Hi replaces Hi, with the small parameter A in- 
troduced to characterize the order of perturbation ex- 
pansion. For the single component interaction Hi{t) — 
S{t)B{t), Eq. ([3]) can be considerably simpUfied. 

One Dimensional Dynamics and the Principle of 
Control. — A primary example is the dynamics, control 
and protection of one normalized state \4>), in the interac- 
tion representation, within the d— dimensional subspace 
Si. [Spontaneous emission, for example, is a case where 
\4)) is an energy eigenstate]. In general |0) is a superpo- 
sition of eigenstates, rather than a single eigenstate, and 
we can rearrange the bases of Si so that \(f>) is one of the 
new orthonormal basis elements. 

Suppose that the initial state 7/(0) = |0) {4>\. The sub- 
system evolves according to the closed equation @ with 
d — 1 and, at time t, ri{t) — b[t) \4>) {<j)\, where, in general, 
b{t) is written as 



h{t) =. exp(-L(i)), 



(4) 



with b{t) < 1 or L{t) > 0. Substituting ri{t) into the mas- 
ter equation ([3]) with d = I, gives an analytic expression 
second order in A, 



L{t) = y / dsC{s), 
Jo 



(5) 



where 



C{s)^ f ds'Y,[&Ms,s-s')<i>^p{s') + h.c.], (6) 

•^0 a/3 



ASc,{s)ASf,is') . Here ASa,{t) 
I !</'). *a/3(i 



and &a/3{s, s') 

Sa{t) - Sa{t) and ^Jt) = 
s) —tTB[Ba{t)Bp{s)pB] is the bath correlation function 
for multi-term system-bath interactions. Note that C(s) 
is a linear function of matrix elements <&q/3. We term the 
time-dependent L{t) a leakage function, by analogy with 
the decoherence function [15[. It describes the leakage 
from Si due to the bath and into the space outside of Si . 
Higher than second order effects in the leakage function 
are included in Eq.([2]) 

L{t) is a functional of the initial state \<p) and any 
added control He, the latter generally through incident 
external fields. Given a time t, the solution of the varia- 
tional equation 6L{t) — with respect to the state in 
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the absence of He yields, to second order, a self protected 
state. Alternatively, solution to this variational equa- 
tion with respect to the incident electromagnetic fields 
for fixed {(p) provides optimal control fields, to second 
order, to protect |(/)) against decoherence. Later in this 
letter we address optimizations with respect to the con- 
trol fields for realistic molecular systems. Optimizations 
with respect to |^) are currently under study. 

As in all approximation techniques, the utility of the 
second order approximation [Eq. ([3])] is examined by 
comparison with exact cases. We consider two examples. 

Example I: Pure leakage. Consider a pure leakage 
case, where the system is a one-dimensional Harmonic 
oscillator in which there is no system-bath interaction. 
The system is described by Hs = cua^a and is polar- 



ized by the interaction Hj = A(i 



Ib, where Ib 



is the unit operator of a bath and a^(a) is a bosonic 
creation (annihilation) operator for the harmonic oscilla- 
tor. For the case of the ground state \(p) = |0), the ex- 
act solution is b^At) = (0| psit) |0) exp(-^ sin^ ^). 

O and dS])] gives the 
, the exact analytical 



The second-order solution [Eqs. 
same result. When 16) ~ \1) 



solution 6cx(i) - 
while b{t) = exp(— 
case of !</.) = ^(|0) 



1 - 

12A^ 



4y 



sm 

,2 ujt 



2f )2exp(-^sin^f ), 
sin^ ^). In the superposition 



bit) 



exp(-i^(sin2f 



-,4 a^t^ 



sm' ^)) and the exact solution 6cx(i) = [1 ^ l^i^ 



A^)sin4f ]exp(- 



4A- 



sin^ ^ ) . Their second-order expan- 



sions are all the same and are essentially equal to one 
another for ^ < 0.1. 

Example II: Spin-bath model within the rotating wave 
approximation. Here the system Hamiltonian reads 



Hq = ea^ + uia^a and Hj 
cr± = (cr^ ± iay)/2. Let |(^) 



a a'), where 



= |0) and let the bath be 
in the state a'f |y) , where |0) is the spin-down state and 
\V) is the vacuum bath state. Solving the problem ex- 
actly gives bex{t) = 1 — sin^(iy^(^ 
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)2-|-A2)- 



The second-order solution to the master equation is 
b{t) = exp(-( 



T sm 



2 (e-cj)t ^ 



-2 — -), which agrees with the 
exact result in second order. They can also be shown 

^ <0.1. 



numerically to be in good agreement when (^g_i^yi 

Quantum Control: Within the framework outlined 
here, the goal of quantum control becomes: given a time 
t of interest, we computationally seek the solution of the 
variational equation dL(t) — with respect to He, with 
the inclusion of any physical constraints of the control. 
The result is the control functional Hc{s'). Note that this 
approach does not just optimize C(s) itself at time s, but 
rather includes the history of the time evolution of C(s). 

Although idealized BB control provides a possible 
mathematical solution, of L(t) = 0, it requires unrealistic 
(zero-width) pulses. Hence, our focus is to replace ideal- 
ized control by an approximate variational or numerical 
solution that minimizes L{t) under realistic pulse energy 



and pulse width constraints. 

Harmonic Oscillator coupled to a harmonic bath.- As 
an application of this framework, consider as the system 
the harmoni c ap proximation (frequency fl) to a Morse os- 
cillator 3, l3| for, e.g., molecular iodine, which allows 
us to employ a simplifying symmetry. The total Hamil- 
tonian, where the bath has £ oscillators, is H = Oa^a + 
t„ . a. " -J2 + a])- The 



-,-int. 



*^^*at and 



X]j=i + XSB, where B 

interaction is separable with S{t) = 
the bath correlation function is $(i) = + 
n(wj)e"*"^* -\- n(a;j)e''^^*} where nj = l/[exp(/3wj) - 1]. 
At low energy, I2 vibrational motion is harmonic, with 

n = 213.7cm-i and a{ujj) = \/ "'" "ioV/ "' " ^it^i ^3 = 
— ^ln(l — '^^^^l — -) , where cod is the cut-off frequency 

at j = e. 

As an example, consider a superposition of the eigen- 
states of the system harmonic oscillator as the state in 
need of protection in a bath at room temperature. Such 
states would be of interest, for example, in pump-dump 
coherent control scenarios [JQ where this is the initially 
pumped state. In that case one would be interested in 
maintaining this state over time scales of ~ 700 fs, the 
system decoherence time[l9|. 

Figure lA shows L{t) for the parameters shown in the 
figure caption. In some cases, S(s, s — s') may be a func- 
tion of (s — s') only, in which case L{t) cx t, shown as 
the dashed line in Fig. lA. This linear dependence on t 
is similar to that of the usual decoherence function [l^ 
in the long time limit. However, at shorter times the 
dependence of 6(s, s — s') on s is seen to contribute non- 
negligibly, leading to oscillatory L{t). Below we will nu- 
merically study the behavior of the leakage function in 
the short-time region to determine the extent to which 
leakage can be controlled. 

Control.- Physically, the origin of the control is that the 
frequency of the system (here, an harmonic oscillator) is 
periodically, dynamically. Stark shifted by the alternat- 
ing field. To this end we employ the realistic control 
Hamiltonian Hs(t) = {n -\- f{t))a^a, which results from 
strong laser pulses acting on electrons that induce an 
additional time-dependent nuclear potential. We model 
the control function as a periodic rectangular interac- 
tion: f{t) =0 for regions other than nr — A < t < nr, 
n integer. Inside these regions f{t) is defined so that 
00 = /r-A •^(^)'^^' That is, for nonzero A, f{t) = 0o/A 
over the control interval, and for A = 0, f{t) — 5{t~nT). 
The functional form contains three main control param- 
eters: the time interval r, the pulse width A and the 
interaction intensity 0o- For comparison with realistic 
pulses, we show L[t) with ideal impulsive phase modu- 
lation (A = 0) in the three lower curves in Fig. 1 (A). 
Clearly, the shorter the control interval, the better the 
control. 

Figure 1(B) shows L(t) with fixed r and (^Oj but with 
different pulse widths A. The results show that the 



4 




FIG. 1: (A) (upper solid curve) L{t) in units of for an 
Harmonic oscillator, in the initial state, 10) = "^(|0) + 
coupled to a bath of 1000 harmonic oscillators at T =300K. 
Here tOd = 0.01(/s)^^ ~ 1.50. The dashed curve shows the 
dominant part, as discussed in the text. The three lower solid 
curves show L(t) with A = 0, 0o = and r = r^^, 
and in the order from bottom to top. (B): L{i) with 
r = and <j)o — n, but different values of A. Solid curve 
A = t/5; Dashed curve: A — r/2; Dotted curve: A — t. (C): 
L{t) with T = and A — ^ but different pulse intensities, 
(j>o — Yfj'f'f from top to bottom. (D): two lower curves, 
L{t) with 00 = tt/IOO, r = ro/100 {tq = ^) and 00 = tt/IO 
and r — ro/10 all with A — t/2. They have the same value of 
fie- The upper curves have the same values of 0o and t but 
different signs of 0o. 



quality of the control is only weakly dependent on the 
pulse width. For example, the control is excellent even 
if the width of the pulse is equal to the control inter- 
val T. In this case the control is equivalent to adding 
a constant frequency fij. = (j)o/T to the harmonic oscil- 
lator frequency, i.e. shifting the system frequency by 
A/r = 1 means shifting the system frequency to (0-1- 
fic)- If (O-t-f^c) > (the cut-off frequency of the bath), 
which is the case in this figure, the function 6(s, s — s') 
oscillates faster than the rate of decay of the bath ^(s). 
The integral, C(s), of the product of the two functions 
generally oscillates around zero so that L{t) is reduced at 
any time t. The way to achieve this goal (e.g., see Ref. 
[l3|) is to increase the interaction of the pulse 4>o or de- 
crease the pulse width t in order to increase flc = 4>q/t. 

The dependence on the intensity is also of interest, as 
shown in Fig. 1(C). The quality of control is seen to 
decrease with decreasing intensity. 

Finally, we consider the effects due to the different 
signs of 00- The two lower curves in Fig. 1(D) cor- 
respond to results for different 0o (or r) but the same 



f^c = 4'ol'^- The upper curves have the same values but 
different signs. They show similar suppression, imply- 
ing that the quality of suppression depends primarily on 
ric- Hence, the above discussions are also valid for the 
negative values of 0o as shown in the two upper curves. 

For the case of a diatomic molecule, we note that the 
AC Stark effect induced by an external laser field inter- 
acting with the electrons decreases (or increases) $7, an 
effect termed "bond softening (or hardening)". Ah initio 
calculations [l^l show that, in the softening case, the fre- 
quency can be reduced by ten percent for in a strong 
laser field. In the case of hardening, if the frequency can 
be enhanced by p percent so that \i {1 -\- p)U, > LOd, then 
leakage control will be effective. This is expected to be a 
considerable technical challenge. 

Conclusion.- We have utilized an exact trace- 
nonpreserving master equation for dynamics of a sys- 
tem subspace, and introduced a leakage function that de- 
scribes leakage from a subspace of interest. We used the 
second-order equation to analyze the quantum dynamics 
of the leakage function, especially for the generic case of 
a harmonic oscillator coupled to a bath of harmonic oscil- 
lators. The realistic pulses required to suppress leakage 
in the presence of the bath are given by the well studied 
finite- width pulses causing AC Stark shifts [H, . A re- 
markable result is that, in the one-dimensional case, one 
can define a leakage function that provides a complete 
description of quantum storage of a general superposi- 
tion state 10). Since quantum storage is of great interest 
at present, the leakage function provides a most conve- 
nient tool for studying the control of quantum states in 
the short time regime relevant to quantum information 
processing and to coherent control. 
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